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ABSTRACT 


Aims: The true nature of primordial magnetic fields (PMFs) and their role in the formation of galaxies still remains elusive. To shed 
light on these unknowns, we investigate their impact by varying two sets of properties: (i) accounting for the effect of PMFs on 
the initial matter power spectrum, and (ii) accounting for their magneto-hydrodynamical effects on the formation of galaxies. By 
comparing both we can determine the dominant agent in shaping galaxy evolution. 

Methods: We use the magneto-hydrodynamics code RAMSES, to generate multiple new zoom-in simulations for eight different host 
halos of dwarf galaxies across a wide luminosity range of 10? — 10° Lo. These halos are selected from a ACDM cosmological box, 
tracking their evolution down to redshift z = 0. We explore a variety of primordial magnetic field (comoving) strengths B, ranging 
from 0.05 to 0.50 nG. 

Results: We find magnetic fields in the interstellar medium not only modify star formation in dwarf spheroidal galaxies but also 
completely prevent the formation of stars in less compact ultra-faints with halo mass and stellar mass below ~ 2.5 - 10? and 3 - 
10° Mo, respectively. At high redshifts, the impact of PMFs on host halos of dwarf galaxies through the modification of the matter 
power spectrum is more dominant than the influence of magneto-hydrodynamics in shaping their gaseous structure. Through the 
amplification of small perturbations ranging in mass from 107 to 10? Mo in the ACDM+PMFs matter power spectrum, primordial 
fields expedite the formation of the first dark matter halos, leading to an earlier onset and a higher star formation rate at redshifts 
z » 12. We investigate the evolution of various energy components and demonstrate that magnetic fields with an initial strength of 
B, 2 0.05 nG exhibit a strong growth of magnetic energy, accompanied by a saturation phase, that starts quickly after the growth 
phase. These trends persist consistently, regardless of the initial conditions, whether it is the classical ACDM or modified by PMFs. 
Lastly, we investigate the impact of PMFs on the present-time observable properties of dwarf galaxies, namely, the half light radius, 
V-band luminosity, mean metallicity and velocity dispersion profile. We find that PMFs with moderate strengths of B; x 0.10nG 
show great agreement with the scaling relations of the observed Local group dwarfs. However, stronger fields lead to large sizes and 
high velocity dispersion. 


Key words. Primordial magnetic fields, Magnetic Fields, Cosmology, Dwarf galaxy, Ultra faint dwarf galaxies, Galaxy evolution, 


Star Formation, Magnetohydrodynamical simulations 


1. Introduction 


Observations reveal the pervasive nature of magnetic fields in 
our Universe, present at all cosmic scales probed so far; be it 
the small scale of planets, stars and galaxies (Stevenson 2010; 
Reiners 2012; Beck 2001; Beck & Wielebinski 2013), or at 
the large scale of galaxy clusters (Clarke et al. 2001; Govoni 
& Feretti 2004; Vogt & EnDlin 2005). In particular for galax- 
ies, magnetic fields play an important role in the process of 
star formation (Padoan & Nordlund 2011; Zamora-Avilés et al. 
2018; Krumholz & Federrath 2019). They have the potential to 
increase the gas fragmentation (Inoue & Yoshida 2019), slow 
down the circumgalactic gas flows (van de Voort et al. 2021), 
suppress the recycling outflows (Grgnnow et al. 2018), and 
shrink the size of galaxies (Martin-Alvarez et al. 2020). The 
role of magnetic fields in shaping the complex multi-phase inter- 
stellar medium is widely recognized (Iffrig & Hennebelle 2017; 
Kórtgen et al. 2019). Observational data from neighboring galax- 
ies indicate that magnetic energy permeates all phases of the 
interstellar medium and is approximately in equipartition with 
the thermal and kinetic energy components (Beck 2015b; Lopez- 
Rodriguez et al. 2023). 


Various theories have been proposed to explain the existence 
of strong galactic fields with wG amplitude that are coherent on 
kpc scales and are detected through radio spectropolarimetry in 
spiral and ultraluminous infrared galaxies (Fletcher et al. 2011; 
Beck 20152; Robishaw et al. 2008; McBride & Heiles 2013), and 
are even present in high-redshift galaxies (Bernet et al. 2008; 
Mao et al. 2017; Geach et al. 2023). Amongst them, dynamo 
mechanisms in collapsed objects (Ichiki & Takahashi 2006; Ryu 
et al. 2008; Naoz & Narayan 2013; Schober et al. 2013; Martin- 
Alvarez et al. 2022), as well as magnetized flows in supernovae 
explosions (Widrow 2002; Hanayama et al. 2005; Safarzadeh 
2018) are only a channel of amplification and require a seed 
magnetic field. During the process of galaxy formation, the col- 
lapse of gas can result in the entanglement of magnetic flux, 
which may also serve as another means of amplification and lead 
to galactic fields (Sur et al. 2012). However, in this case, a strong 
micro gauss initial seed, comparable to what is observed today, 
is required. 


Seed magnetic fields for the galactic dynamo can be gener- 
ated by specific astrophysical mechanisms, like the Biermann 
battery during structure formation (e.g. Kulsrud et al. 1997; 
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Ryu et al. 1998) or reionization (e.g Subramanian et al. 1994; 
Gnedin et al. 2000) or the Weibel instability (Sironi et al. 2023; 
Zhou et al. 2024). Alternatively, cosmic magnetic fields can be 
remnants from the early Universe, generated in the course of 
phase transitions (Hogan 1983; Ratra 1992a; Ellis et al. 2019) 
or during the inflationary expansion (Turner & Widrow 19882; 
Adshead et al. 2016; Domcke et al. 2019; Fujita & Durrer 
2019). Inflation, in particular, provides an ideal mechanism for 
producing large-scale magnetic fields (Ratra 1992b; Turner & 
Widrow 1988b). These fields could undergo further amplifica- 
tion through a small-scale chiral dynamo (Joyce & Shaposh- 
nikov 1997; Schober et al. 2022). Such primordial magnetic 
fields, when generated during this early epoch, would be am- 
plified during the collapse of density perturbations, and directly 
explain the magnetic fields observed in galaxies (Kandus et al. 
2011; Martin-Alvarez et al. 2021). 


Observational evidence of non-negligible magnetism even in 
the intergalactic medium voids has recently sparked interest in 
this intriguing possibility (Neronov & Vovk 2010; Taylor et al. 
2011; Tavecchio et al. 2011). These observations are based on a 
missing y-ray signal from TeV blazars located behind the voids. 
The presence of magnetic fields in these relatively quiet environ- 
ments is difficult to explain purely by an astrophysical process in 
the late Universe (Bertone et al. 2006; Furlanetto & Loeb 2001; 
Vazza et al. 2017; O'Sullivan et al. 2020; Arámburo-García et al. 
2021; Tjemsland et al. 2023), and would perhaps favor a pri- 
mordial origin. The recent discovery of magnetized filaments 
between galaxy clusters reinforces this scenario (Govoni et al. 
2019; Carretti et al. 2022). For comprehensive reviews on differ- 
ent scenarios for the generation of magnetic fields see Widrow 
(2002), Kandus et al. (2011), and Subramanian (2016). 


The main interest in magneto-genesis scenarios with a pri- 
mordial origin stems from their capacity in simultaneously ex- 
plaining the = 10716 G intergalactic magnetic fields as well as 
their ~ 10^? G galactic counter-parts without imposing a time 
constraint for rapid amplification and subsequent reorganization 
into their observed small and large-scale ISM distribution (e.g., 
Beck 2015b; Borlaff et al. 2023). The existence of such mag- 
netic fields in the pre-recombination plasma may even help to 
alleviate the Hubble tension (Jedamzik & Pogosian 2020). This 
result emerges from an increased rate of recombination in the 
presence of strong primordial magnetic fields (Galli et al. 2022). 
Before and during the Epoch of Recombination, these primor- 
dial magnetic fields generate density perturbations in addition to 
those produced by inflation. These perturbations are mapped on 
the last scattering surface of the cosmic microwave background 
(CMB) photons, where primordial magnetic fields can imprint 
a variety of signals. The trace of a stochastic primordial mag- 
netic field on the CMB map could be observed in the excess 
temperature and polarization anisotropies (Subramanian et al. 
2006; Durrer 2007; Widrow et al. 2012; Planck Collaboration 
et al. 2016b). These magnetically-produced anisotropies display 
a strongly non-Gaussian pattern (Caprini et al. 2009; Seshadri & 
Subramanian 2009; Trivedi et al. 2014). A comprehensive anal- 
ysis of the CMB constrains the primordial magnetic field to have 
an upper limit for its amplitude of a few nano Gauss, favour- 
ing a nearly scale invariant spectrum smoothed over the spatial 
scale of 1 Mpc. Alternative observational constraints of primor- 
dial magnetic fields stem from Big Bang nucleosynthesis (Ker- 
nan et al. 1996, which provides the strongest constraints on all 
length scales), the Sunyaev-Zel'dovich statistics (Tashiro et al. 
2009), the propagation of ultra-high energy cosmic rays (Bray 
& Scaife 2018; Alves Batista & Saveliev 2019), the two-point 
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shear correlation function from gravitational lensing (Pandey & 
Sethi 2012), and Lyman-a forest clouds (Pandey & Sethi 2013). 

Density perturbations produced by primordial magnetic 
fields influence the formation of first structures and the total mat- 
ter power spectrum. This impact commences when the ionized 
matter is decoupling from photons (Trivedi et al. 2012). At this 
stage, magnetic Alfvén waves can produce additional motions 
in the nearly homogeneous background (Tsagas & Maartens 
2000). These additional density and velocity perturbations in 
the baryonic matter (Kim et al. 1996; Wasserman 1978b) are 
gravitationaly coupled to dark matter over-densities, and con- 
sequently modify the number distribution of the first dark mat- 
ter halos (Varalakshmi & Nigam 2017; Cheera & Nigam 2018). 
The magnetically-produced perturbations dominate the standard 
A cold dark matter (ACDM) power spectrum at large k-modes 
(Wasserman 1978b; Gopal & Sethi 2003; Shaw & Lewis 2012), 
and therefore affect the formation of small scale structures. This 
region of k-space encompasses the low mass dark matter halos, 
hosting the smallest galactic systems, dwarf galaxies (see Simon 
2019, for a recent review). 

Despite being the most abundant type of galaxies, dwarf 
galaxies are faint and difficult to observe. Therefore, the best 
studied population of dwarf galaxies have been discovered in 
our proximity (McConnachie 2012a, updated online catalog), 
and in the nearby Virgo cluster (Binggeli et al. 1985). With the 
new window the James Webb Space Telescope (JWST) opens 
to the young Universe, high redshift dwarf galaxies further ex- 
pand our observational data (Jeon & Bromm 2019; Gelli et al. 
2021; Yajima et al. 2022). The existing observations reveal the 
scaling of dwarf properties, from brighter systems, to the regime 
of ultra faint dwarfs. Some scaling relations that characterize 
the properties of dwarf galaxies are the size-luminosity (e.g. 
Patel et al. 2018), metallicity-luminosity (Lequeux et al. 1979; 
Skillman et al. 1989; Garnett 2002; Tremonti et al. 2004), and 
luminosity-velocity dispersion (e.g. Battaglia et al. 2008; Walker 
et al. 2009; Fabrizio et al. 2011, 2016). Thus far cosmological 
simulations have been relatively successful in reproducing the 
scaling relations for dwarf spheroidals (Jeon et al. 2017; Mac- 
ció et al. 2017; Revaz & Jablonka 2018a; Escala et al. 2018; 
Wheeler et al. 2019; Agertz et al. 2020; Applebaum et al. 2020; 
Prgomet et al. 2022). However, the assembly of fainter galaxies 
within dark matter halos is still subject to persistent tensions in a 
ACDM paradigm. The most prominent are the diversity of sizes 
(Oman et al. 2015; Relatores et al. 2019; Revaz 2023) and diver- 
sity of rotation curves (Oman et al. 2015), the satellite planes 
(Pawlowski & Kroupa 2013), and the too-big-to-fail problem 
(Boylan-Kolchin et al. 2011, 2012). For a recent review on the 
long-standing tensions between ACDM and Local Group obser- 
vations see Bullock & Boylan-Kolchin (2017) and Sales et al. 
(2022). 

Modern cosmological simulations have made significant 
progress in addressing small-scale challenges in the standard 
model by incorporating more precise treatments of baryonic 
physics (e.g., Gutcke et al. 2021; Koudmani et al. 2022; Martin- 
Alvarez et al. 2023; Hopkins et al. 2023). The importance of 
magnetic fields in shaping galaxies and their the interstellar 
medium has been recognized and studied by a variety of simula- 
tions (Martin-Alvarez et al. 2020; Pakmor et al. 2017; Su et al. 
2020). However, accounting for additional factors such as pri- 
mordial magnetic fields, which are often ignored in models, will 
provide new insights into existing tensions. Presenting a suite of 
cosmological simulations, Sanati et al. (2020) showed the impact 
of primordial magnetic fields on the formation of dwarf galaxies 
during the Epoch of Reionization. The aim of this manuscript is 
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to expand on that work, by investigating the role played by pri- 
mordial magnetic fields in the formation of dwarf galaxies and 
shaping their global properties, conducting a set of cosmolog- 
ical zoom-in simulations that are designed to self-consistently 
model: (1) the formation and evolution of dwarf galaxies emerg- 
ing from a magnetically-distorted ACDM matter power spec- 
trum, (ii) the growth and amplification of uniform primordial 
magnetic fields. 

This paper is organized as follows. Section 2 reviews the 
role played by primordial magnetic fields in the early history 
of the Universe, and their influence on the ACDM matter power 
spectrum. The numerical framework to generate and evolve our 
simulations is described in details in Section 3. The results are 
presented in Section 4. We first discuss the first stages in the 
evolution of dwarf galaxies, including the first star forming gas 
halos in Section 4.1. Then in Section 4.2 we present the evolu- 
tion of different energy components. In Section 4.3 we follow the 
star formation history and the scaling relations of our simulated 
dwarf models. Finally, a summary of our main conclusions and 
a brief discussion is presented in Section 5. 


2. Inclusion of primordial magnetic fields in ACDM 


In Section 1, it was discussed that cosmic magnetic fields may 
have originated from the inflationary epoch or phase transitions 
in the early Universe. Such primordial fields with a stochastic 
nature modify the baryonic clumping factor, thereby affect the 
formation of initial structures and the ACDM matter power spec- 
trum. This section provides a detailed review of such effects. 


2.1. Primordial magnetic field power spectrum 


Primordial magneto-genesis scenarios motivated by inflation, 
conventionally result in a homogeneous, isotropic and Gaus- 
sian random initial magnetic field B (x,t) (Wasserman 1978a; 
Kim et al. 1996; Gopal & Sethi 2003). The two-point correla- 
tion function of a non-helical magnetic field B in Fourier space, 
therefore, can be described by the relation (Landau & Lifshitz 
1987): 


rut ) 


(Bi(k) Bi(k’)) = (27) à(k — k') ——— Pa), (1) 


where k = |k| is the comoving wave number. The power spec- 
trum of the primordial magnetic field Pg(k) is defined by a sim- 
ple power law: 


Pp(k) = Agk”, (2) 


and projected into the Fourier space, orthogonal to k;, through 
the transverse projection tensor! P; i(k) = ójj — kik; /k2. Here ng 
is the slope, and Ag the amplitude of the magnetic power spec- 
trum. Commonly, the amplitude of the spectrum is defined as the 
variance of magnetic field strength at the present time B; (x, to), 
smoothed at scales of A = 27/k, = 1 Mpc, as 


ru ) 


(Gs) = -3 f dk k? 2x 6(k — k') ——— Pg(k), (3) 


where the integral upper limit ke = ae From Eq. 3 the 
amplitude Ag is expressed as (Shaw & Lewis 2012) 
(QQays*? B 


m — (4) 
2p (22) Ke 


! The importance of the projection tensor lies in the fact that estimating 
the perpendicularly projection to the line of sight is easier than measur- 
ing the tensor field from the observations in large-scale structure. 


With this definition, ng and B; fully characterize the primordial 
magnetic fields. They constitute our main free parameters. In the 
following section, we will describe the incorporation of primor- 
dial magnetic fields into the evolution of density perturbations 
using the linearized Newtonian theory. 


2.2. Evolution of density fluctuations 


In the ACDM framework, all cosmic structures are the result of 
primordial quantum field fluctuations. In the early Universe, the 
causal connection between these small over-densities shapes a 
nearly scale invariant power spectrum P(k) o k”, as a function 
of wave number k, with n « 1 (Harrison 1970; Zeldovich 1972; 
Peebles & Yu 1970; Planck Collaboration et al. 2016b). 

In the linear Newtonian paradigm, the evolution of ini- 
tial density perturbations in the baryonic matter 6,(x,f) = 
ópp(x,t)/pp and the collisionless dark matter ópw(x,f) = 
dppm(*, t)/ppw is described by the two following coupled equa- 
tions (see, e.g., Dodelson 2003), 


O76, 4p; 0ó, 
+12 A(t z Ne E 
Or me 3py | at 
av 5p + 47G [ppdp + ppwópw] i 
and 
026; 
ae = 4nG [ps6 + PpMOpM | (6) 


Here a is the scale EN H = a/a is the Hubble rate, pp and 
py are respectively the baryon and photon mass density, n, is the 
electron number density, or is the Thomson cross section for 
electron-photon scattering, and c; is the baryon sound speed. 

Prior to recombination, due to the high rate of photon-baryon 
scatterings, radiation and baryonic matter behave as a sole tightly 
coupled conducting medium, with the magnetic pressure be- 
ing insignificant compared to the pressure of the photon-baryon 
fluid. As soon as photons start to free stream the evolution of 
baryonic perturbations can no longer be adequately described by 
Eq. 5. It is thus necessary to incorporate compressional and rota- 
tional perturbations in this ionized medium. For that the Lorentz 
force can be formulated using the tangled field, 
segs. ee (7) 

Anpp(to) 

We can now solve the magnetohydrodynamics equation de- 
scribing the growth of total matter density perturbations (see, 
e.g., Subramanian & Barrow 1998; Sethi & Subramanian 2005), 


o, S(x,t 
= = 42Gpndm + = SOM (8) 
E Pm @ 
where óm = (QpMÓpM +Pbôb)/Pm, With Pm = (Cpm+pp) stands 


for the summation over both baryonic and dark matter perturbed 
components. For scales above the Jeans mass, the term represent- 
ing fluid pressure, cV óp, can be neglected. Following recombi- 
nation, the damping term in Eq. 5 attributed to the radiative vis- 
cousity can also be neglected, as the damping effect caused by 
the Hubble expansion becomes more important (Jedamzik et al. 
1998; Subramanian & Barrow 1998). The growth of perturba- 
tions generated by any inhomogeneous? magnetic field is then 
found to be 


2 The homogeneous solutions are related to fluctuations produced by 
sources before recombination. 
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Fig. 1: The impact of primordial magnetic fields on the matter 
power spectrum varying with increasing field strength. Dashed 
lines illustrate the distribution of matter induced by primordial 
magnetic fields, while the solid lines represent the combination 
of these magnetic-induced matter perturbations and the perturba- 
tions produced during inflation. The unperturbed ACDM spec- 
trum is shown as the dark blue solid line. The amplitude of the 
magnetic field is varying from 0.05 to 0.50 nG, while the spec- 
tral slope is kept constant at ng = —2.9. The triangles with error 
bars represent Lyman Alpha data (Tegmark & Zaldarriaga 2002) 
as the observational constraints on the matter power spectrum at 
the smallest mass scales. 


3t rec 3 


rd S(x, trec 
rec (x ) Pb zis 3) (9) 


Om,PMFs(X, t) = 
aÈ (trec) Pm 


9 [ 1 VP f 
10 \ frec 5t 2 

where frec is the time of recombination epoch. The spatial 
variation of ôm pmrs(x, t) in Eq. 9 can be tracked through the mag- 
netic source term S (x, t). Therefore, it is expected that the power 
spectrum of the magnetic field directly affects over-densities on 
various scales, and ultimately influences the total matter power 
spectrum. 


2.3. Impact on the total matter power spectrum 


To understand the effect of primordial magnetic fields on the 
ACDM matter power spectrum, computing the ensemble aver- 
age of density perturbations is required. In a primordially mag- 
netized early Universe, the total perturbations are the sum of 
those generated by inflation and the additional perturbations in- 
duced by primordial magnetic fields (Kim et al. 1996; Gopal & 
Sethi 2003). With a Fourier space expression, it leads to 


P(k, t) = (lôn (k, £) + ôm ps. D] Onk, D) + On ps (e, D]. (10) 


Figure 1 showcases P(k,t) at the present epoch. The dark 
blue line demonstrates the ACDM matter power spectrum, re- 
sulting from the ensemble average of [6,,(k, t)|?. This represents 
density perturbations obtained from the coupled equations 5 and 
6 in the absence of primordial magnetic fields. To obtain the 
matter power spectrum induced by primordial magnetic fields, 
which is demonstrated by dashed lines, we need the Fourier 
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transform of Eq. 7. Because of the two spatial derivatives in 
S (x,t), this results in Ppmps(k, t) = ([ón purs (Kk. 11?) oc K* BZ ox 
Ap k?"»*". At large k modes, perturbations produced by a scale- 
invariant magnetic spectrum (i.e. ng = —3 and Ppmps(k) ~ k) 
dominate the inflationary induced perturbations. However, the 
magnetic pressure leads to an increase of the Jeans length and 
denotes the lower limit at which perturbations can collapse. Kim 
et al. (1996) suggested a formulation for the magnetic Jeans scale 
which occurs at the transition between growing modes and os- 
cillatory modes, 


| 570p VG 


B, uu 


kjg 
By considering the magnetic Jeans length, we arrive at an expres- 
sion for the shape of dashed lines in Fig. 1, which corresponds 
to the matter power spectrum induced by primordial magnetic 
fields, ([ó Pmr (k, 0]?), (Gopal & Sethi 2003), 


Ppmps(k) ~A ponet + Bi"? + C K+! k6 qus 


max max 


(12) 


Here A, B and C are coefficients which, similarly to Ag, only 
depend on the spectral index. The magnetic counterpart of the 
Jeans length kmax = kjg is regarded as the minimum length scale 
above which the growth of density perturbations is suppressed. 
Equation 12 highlights that for ng « —1.5 perturbations on small 
scales are amplified until they reach kmax, beyond which they are 
rapidly quenched. Finally, solid colored lines in Fig. 1 indicate 
the total power spectrum obtained from Eq. 10 with the inclusion 
of both inflationary and magnetically generated perturbations. 

To generate the power spectra presented in Fig. 1, we use 
a modified version of the CAMB code (Shaw & Lewis 2012) 
that incorporates primordial magnetic fields. This version ex- 
plicitly calculates the non-linear effects of magnetic pressure 
and the non-negligible viscous damping prior to recombination. 
Therefore, there is no need for an artificial cutoff wave num- 
ber, making the power spectrum more accurate compared to an- 
alytical studies in the literature. The range of parameters ex- 
plored in Fig. 1 is carefully selected to affect the power spectrum 
while still adhering to the constraints presented in Section 1. 
Therefore, we use a nearly scale-invariant power spectrum with 
ng = —2.9, and examine a variety of magnetic field amplitudes, 
Bı = 0.05, 0.10, 0.20, 0.50 nG. On the top axis, we show the 
link between the scales in Fourier space and the corresponding 
mass scales in real space. For that, we define the mass enclosed 
within a comoving Lagrangian sphere of radius r; = 27/k at the 
present epoch, (see e.g., Bullock & Boylan-Kolchin 2017), 


HU = OH 3 
3 [Fm 2G 1 


2 3 
Q h ri 
1.71 x 10!! MS | — 
DUE (Ss) ss] (ats) 


where Ho is the Hubble constant. With this definition, it be- 
comes evident that the tilt in the matter power spectrum appears 
in the mass ranges that encompass the dark matter halos of dwarf 
galaxies. Consequently, we expect that the formation and charac- 
teristics of dwarfs could be substantially affected by the presence 
of primordial magnetic fields. 


M, = 


(13) 


3. Numerical methods and simulations 


We generate all the cosmological simulations studied in this 
work employing our modified version of the RAMSES code 
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(Teyssier 2002). In addition to the collisionless dark matter and 
stellar particles, RAMSES employs an adaptive mesh refinement 
octree grid to solve the evolution of gas. The code models ideal 
magneto-hydrodynamics (MHD; Fromang et al. 2006; Teyssier 
et al. 2006) in addition to treating baryonic physics, such as 
redshift-evolving and uniform UV heating, gas cooling, star for- 
mation, and stellar feedback. Below we provide a brief summary 
of its essential elements. 


3.1. Numerical setup 


Ideal MHD In its ideal MHD implementation, RAMSES employs 
a Constrained Transport (CT) method (Teyssier et al. 2006; Fro- 
mang et al. 2006) to solve the equations that govern the evo- 
lution of magnetic fields. The induction equation is solved in a 
conserved integral form on the cell faces. This requires magnetic 
fields to be stored as six fields on the cell faces. This is unlike 
all the hydrodynamic quantities in the simulation, namely, den- 
sities, velocities, and energy components which are stored at the 
center of each gas cell. The CT method ensures that the magnetic 
field has zero divergence up to numerical precision, preventing 
any unwanted modifications of conserved quantities (Tóth 2000) 
or the emergence of magneto-hydrodynamical artifacts (Hopkins 
& Raives 2016). 


In the ideal MHD setup, assuming a highly conductive 
medium, the induction equation is solved with negligible dif- 
fusivity but we note that numerical diffusivity is present because 
finite size of the grid cells. Additionally, the distortion of the 
primordial field due to the velocity perturbations is only signifi- 
cant at sufficiently small scales and can be disregarded at galac- 
tic scales (Wasserman 1978b; Peebles 1980). In the absence of 
non-ideal magnetic sources such as the Biermann battery, this 
leads to ô (a2 B,)/Ót = 0 for the time evolution of the field. As 
the growth of compressional modes is suppressed before recom- 
bination, it solves as By(x,t) = By(X, trec) a7(trec)/a2(t), where 
B(x, trec) refers to the value of the magnetic field at recombina- 
tion (f = frec). In this work, the primordial magnetic field is mod- 
eled by an ab-initio B, that is seeded uniformly and aligned with 
the z-axis of the computational domain. This choice of a uniform 
field is common practice amongst MHD simulations (e.g., Vazza 
et al. 2014; Pakmor et al. 2017; Martin-Alvarez et al. 2018), but 
we note that smaller-scale perturbations for such strong initial 
fields may influence the magnetic field amplification and galaxy 
formation during the early stages of simulations. 


Radiative cooling and heating processes In addition to pri- 
mordial gas cooling, we account for metal-line cooling accord- 
ing to the gas metallicity. Above temperatures of 10^ K we in- 
terpolate the corresponding cLoupy tables (Ferland et al. 1998). 
Below 10^ K we follow fine structure metal cooling rates from 
Rosen & Bregman (1995). We model the process of reioniza- 
tion, based on the prediction from Haardt & Madau (1996), us- 
ing a redshift-dependent UV background which we initialize at 
z = 10. Hydrogen self-shielding against the ionizing radiation 
is incorporated by suppressing the UV-background heating for 
gas densities above ny = 0.01 g cm. To account for the pre- 
enrichment by the Population III stars we assume a metallicity 
floor of [Fe/H] = —5. Above this metallicity, fine-structure line 
cooling of atomic carbon and oxygen (Bromm & Loeb 2003), 
lead to gas fragmentation and the formation of low-mass stars. 


Star formation We model the process of star formation employ- 
ing a magneto-thermo-turbulent (MTT) star formation prescrip- 
tion, presented in more detail in the hydrodynamical version by 
Kimm et al. (2017) and Trebitsch et al. (2017), and in its MHD 
version by Martin-Alvarez et al. (2020). As a brief summary of 
this star formation model, gas is allowed to transform into stel- 
lar particles only in cells that are at the highest level of refine- 
ment, where the local combination of magnetic, thermal and tur- 
bulent support is overcomed by the gravitational pull (Rasera & 
Teyssier 2006). The conversion of gas into stars in star forming 
cells follows the Schmidt law (Schmidt 1959), 
Pgas 


D, = EF . 
Px fg 


(14) 


Here the star-forming gas with a density of pga, is converted into 
stars, with free-fall time £g. In our formulation, the efficiency of 
this conversion eg is not a constant, but rather computed locally 
based on the gas properties of each region. It is calculated fol- 
lowing the multi-scale model of (Padoan & Nordlund 2011) as 
presented in Federrath & Klessen (2012). 


Stellar feedback Each stellar particle in our simulations corre- 
sponds to a single stellar population that is characterized by its 
own initial mass function (IMF). The IMF is modeled as a prob- 
ability distribution function following a Kroupa shape (Kroupa 
2001) and normalized over the complete range of masses. This 
allows each stellar particle of mass 294.81 Mo to stochastically 
populate stars within the mass interval of [0.05 — 50] Mg dur- 
ing the initial 50 Myr of its formation. The number of exploding 
supernovae for each particle is then calculated at each time step 
based on the lifetimes of stars it contained. 

We use the mechanical supernovae feedback prescription of 
Kimm & Cen (2014), for stellar particles undergoing a super- 
nova event. In this approach, the momentum injected by super- 
nova explosions is determined by the physical characteristics of 
the gas being swept up, such as its density ny, and metallicity Z, 


psn(E, nu, Z) = 3x kms Mo ES nq f(Z), — (05) 


where the momentum input is decreasing with metallicity as 
f(Z) = max [Z/Zo, 0.01] 9-4. Here Es; is the energy in the 
unit of 10?! erg, and ny the hydrogen number density. This 
method ensures that the feedback from supernovae is accurately 
modeled at all blast wave stages, from the initial free expan- 
sion to the final momentum-conserving snowplow phase. Along 
with the momentum, energy is deposited into the neighboring 
cells. The specific energy of each supernova (SN) has a value 
Esn = Esn/Msn, where Esn = 10°! erg and Msn = 10 Mo. Each 
supernova also returns a fraction of stellar mass back to the ISM. 
We use nsn = 0.213 for fraction of Msn returned as gas mass, 
and meta = 0.075 for the newly synthesized metals. This set of 
values is selected to match the model galaxies in the fiducial set- 
ting with the observed global properties of Local Group dwarfs 
(See Sec. 4.4). Our feedback prescription does not include stellar 
winds, local radiation from young stars, photo-electric heating, 
and cosmic rays, all of which, in addition to supernovae feed- 
back, can potentially reduce the final stellar mass of our galax- 
ies. 


3.2. Initial conditions 


The initial conditions are generated using the MUSIC code (Hahn 
& Abel 2011), and the cosmology of Planck Collaboration et al. 
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Table 1: Parameters varying in the simulation runs. Columns are as follows: 1) Model ID. 2) Comoving strength of the initial primordial magnetic 


field. 3) Matter power spectrum used to generate the initial conditions. 


ID BOO BOS B10 B20 B50 | BOSps Bl10ps B20ps B50ps 
Ba [nG] 0 0.05 0.10 0.20 0.50 0.05 0.10 0.20 0.50 
Power Spectrum ACDM ACDM+PMFs 


(20162) with Qa = 0.685, O,, = 0.315, Q, = 0.0486, and h = 
0.673. All simulations are started at redshift z = 200, ensuring 
that the rms variance of the initial density field, cs, lies between 
0.1 and 0.2 (Knebe et al. 2009; Ofiorbe et al. 2015). Except for a 
small subset, all of our simulations are run until redshift z = 0. 

The twelve halos studied in this work are selected from dif- 
ferent regions in the dark matter only cosmological box of Re- 
vaz & Jablonka (20182). Those regions are identified to form a 
dark matter halo with a mass typical of spheroidal and ultra-faint 
dwarf galaxies at redshift z — 0. We re-simulate those halos, 
including the full treatment of baryons described above. Using 
the zoom-in technique, for each halo we refine a 3D ellipsoid 
of size about 0.85 cMpc across, positioned in the center of the 
cubic simulation box with Lgox = 5.11 cMpc per side. The size 
of the ellipsoid is determined such that it encompasses all par- 
ticles that eventually reside within each target halo by redshift 
z = 0 (Oforbe et al. 2014). In this refined region, we achieve 
dark matter mass resolution of mpy ~ 4 - 10° Mo. We gradually 
degrade the resolution, from level 10 to 6 of MUSIC outside the 
zoom region ?. 

Initially, the domain is discretized with a uniform grid of 
1024? cells. This resolution is preserved within the zoom region, 
while the grid is de-refined to level 6 elsewhere. Throughout the 
course of the simulation, the adaptive refinement criteria come 
into play to effectively resolve dense and Jeans-unstable regions. 
When the total dark matter and gas mass within a grid cell ex- 
ceeds 8 mpm, or when the size of the grid cell surpasses 4 local 
Jeans length, a parent grid cell is split into 8 equal child cells. 
This process follows the octree structure of RAMSES, where the 
size of cell i is determined by the refinement level of cell /; ac- 
cording to Ax; = 1/ 2! Lox. In our simulations, with a maximum 
refinement level of 20, the initial grid undergoes adaptive refine- 
ment to achieve a minimum cell width of approximately 5 pc. 

These initial conditions and the configuration of the simula- 
tions focus on capturing and resolving the formation of galax- 
ies in halos within the perturbed range of the total matter power 
spectrum in Fig. 1. 


3.3. Set of simulations 


We run nine sets of simulations as listed in Table 1. The run 
which has no magnetic fields and employs a ACDM initial con- 
figuration is our reference, hereafter denoted as BOO. In the first 
series of simulations, a uniform magnetic field is introduced, 
aligned with the z-axis of the simulation box. The field has an 
initial comoving strength varied from 0.05 to 0.50nG, which 
later evolves self-consistently throughout the simulation. The 
initial matter power spectrum is set using the ACDM paradigm, 
where no additional perturbation due to primordial magnetic 
fields is included. Models in this set are denoted by B05, B10, 
B20, B50. In the second series of simulations, in addition to the 
initial uniform field of comoving strength ranging from 0.05 to 


? One resolution level / corresponds to N = (25? particles in the full 
cosmological box. The particle mass is thus decreased by a factor of 
eight between two levels. 
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0.50 nG, the effect of primordial fields on the matter power spec- 
trum is taken into account. In these models, denoted as BO5ps, 
B10ps, B20ps, B56ps, the initial conditions are generated from 
a ACDM+PMFs power spectrum. In each model, the amplitude of 
the primordial magnetic field is the same as the uniform field 
in the MHD configuration initiated at redshift z 2 200, and the 
spectral index is set to ng — —2.9 to have a scale-invariant mag- 
netic power spectrum. 

All the simulations are designed to investigate primordial 
magnetic fields with strengths allowed by current cosmologi- 
cal constraints (see Section 1). From Fig. 1 it appears that a 
stronger amplitude than the maximum one used in this work 
(Ba = 0.50nG) would affect the power spectrum at scales of 
k < 10A Mpc^!. However, such a modification is ruled out by 
the observations of the Lyman-a forest (Tegmark & Zaldarriaga 
2002; Pandey & Sethi 2013). The smallest amplitude employed 
(B4 = 0.05 nG) lies at the threshold where primordial magnetic 
fields start to affect the ACDM power spectrum and potentially 
exert an influence on structure formation (Shaw & Lewis 2010; 
Pandey et al. 2015; Sanati et al. 2020). 


3.4. Extraction of the observables, luminosity and metallicity 


Galactic region Global galactic properties, such as gas mass 
and specific energies, are measured inside the galactic region. It 
is defined using the virial radius of the dark matter halo R209 at 
redshift z = 0. The galactic region is the spherical volume cen- 
tered on the position of the galaxy with radius equal to the co- 
moving evolution of ~ 1.5 Roo (z = 0). This method is employed 
to adaptively increase the size of the studied region as the galaxy 
and its hosting dark matter halo grow in size. In order to obtain 
the centre of the galaxy, we tag the system by finding its position 
at redshift z ~ 7. From that point onwards, at each coarse time- 
step of RAMSES we re-compute the position of the galaxy centre. 
To do this, we simply obtain the new center of mass of the most 
central ~ 50 stellar particles in the previous iteration. We use this 
approach to compute various galaxy properties on-the-fly, with 
time resolution equivalent to that of the coarse time-stepping of 
the simulation. 


V-band luminosity The galaxy V-band luminosity (Ly) is ob- 
tained by summing the luminosities of all stellar particles lo- 
cated within 1.8 times the halo’s half-light radius. Their mass 
is converted into luminosity using the stellar population synthe- 
sis model of Vazdekis et al. (1996) computed with the revised 
Kroupa (2001) IMF. Where necessary, the luminosities are inter- 
and extra-polated in age and metallicity using a bivariate spline. 
We ignore dust absorption, as its role only becomes critical for 
intrinsic magnitudes brighter than Mag ~ —20 (Ma et al. 2018). 


Metallicity Due to the limited number of stars forming in dwarf 
galaxies, selecting a representative mean [Fe/H] is a non-trivial 
undertaking. Since the metallicity distribution is sparsely sam- 
pled, calculating the mode (peak) of the metallicity distribution 
function can result in significant uncertainties. In this study, we 
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Table 2: Global properties of the eight halos simulated in our fiducial hydrodynamical model B00, at redshift z = 0. The first column gives the 
halo ID following Revaz & Jablonka (20182). Ly is the V-band luminosity and M, the stellar mass. Moog is the virial mass, i.e., the mass inside 
the virial radius R99. Cros is the line-of-sight velocity dispersion and [Fe/H] is the abundance ratio of iron with respect to hydrogen. 


Halo ID Ly M, M200 Roo npo cros [Fe/H] 
[IO Lo] [10°Mo] [10°Mə] [kpc] [kpc] [km/s] [dex] 
h025 9.55 19.67 8.90 65.6 0.57 9.4 -1.95 
h070 1.02 2.42 2.16 40.2 0.50 6.0 -2.73 
h063 0.58 1.44 2.49 42.2 0.82 4.4 -2.41 
h177 0.31 0.78 0.61 26.4 0.45 4.5 -2.42 
h111 0.16 0.38 1.25 33.6 0.35 5.0 -2.94 
h159 0.24 0.59 0.78 28.7 0.42 4.8 -2.69 
h277 0.09 0.21 0.90 30.0 2.00 5.9 -2.65 
h190 0.06 0.15 0.57 25.8 0.72 3.8 -3.04 


use the method proposed in Sanati et al. (2023) to determine the 
mode based on a fitting analytical formula derived from a simple 
chemical evolution model. The error on [Fe/H] is taken at the 
maximum of the errors obtained by the different methods, peak 
and mode. 


Velocity dispersion The line-of-sight stellar velocity disper- 
sion, Gros, is calculated for seven different lines of sight inside 
a 1 kpc cylindrical radius. The value quoted for each galaxy rep- 
resents the mean of these values. 


4. Results 


Our simulations follow the evolution of six different ultra-faints 
and two spheroidal dwarf galaxy with final luminosities ~ 10? — 
10? Lo. Table 2 summarizes the main properties of model galax- 
ies in our fiducial hydrodynamical model BOO, at redshift z = 0. 
IDs are those of Revaz & Jablonka (2018a, see their Table 1). 
The seven other columns show the galaxy total V-band lumi- 
nosity Ly, total stellar mas M,, virial mass M99, virial radius 
R209, mean stellar velocity dispersion o1 os and peak value of 
the metallicity distribution function [Fe/H]. 


4.1. Impact of magnetic fields on gas distribution 


Figure 2 shows the density-weighted projection of gas density 
in the simulated galaxy h070 (see Table 2). The magnetic field 
strength is increasing from 0.05 nG in the leftmost to 0.50nG in 
the rightmost panel. Top and bottom panels correspond to ACDM 
and ACDM+PMFs simulation setup, respectively. Snapshots are 
taken at redshift z = 9, which corresponds to the intense phase 
of star formation in dwarf galaxies (Okamoto et al. 2012; Brown 
et al. 2014; Sacchi et al. 2021; Gallart et al. 2021). The snapshot 
in model B50ps is shown at redshift z = 13 where it stopped due 
to computational expenses. 

We compute the gas mass in each snapshot inside the galactic 
region (see Section 3.4) which is illustrated by the white circles. 
The galactic region covers the main halo of the targeted dwarf 
galaxy and its surrounding in a 1.5 x Ryjr radius. Mga in the top 
left of each panel, shows the gas mass within the galactic region 
averaged over the first 2 Gyr, before the dwarf is completely de- 
pleted of gas. 

The upper panels correspond to the first set of simulations, 
where the magnetic fields are initialized as uniform values, as 
described in Section 3, and the matter power spectrum is that of 
the classical ACDM paradigm. These simulations are designed 


to explore the contribution of magnetic fields in shaping the in- 
ternal structure of dwarf galaxies. Here, magnetic fields generate 
additional pressure support. Along with other hydrodynamical 
pressures such as the thermal or kinematic, this counteracts the 
gravitational field and effectively hampers the collapse of gas. 
It particularly affects the first stages of galaxy formation, i.e., 
for z > 6. The result is that the gas collapse in the first miniha- 
los is delayed by several hundred Myrs. The gravitational field 
in minihalos is therefore less likely to retain the gas when it is 
heated and pushed out by the UV-radiation at lower redshifts. 
When B; = 0.50nG the collapse of gas is completely halted, as 
shown in the most right panel. This effect leads to a reduction 
trend in the averaged gas mass with increasing B. 


The bottom panels of Fig. 2 represent the second set of simu- 
lations, where in addition to varying the strength of the magnetic 
field seeds, we include their effect on the matter power spec- 
trum. These simulations are designed to explore the net effect 
of primordial magnetic fields through the altered matter power 
spectrum on the one hand, and the direct magnetically-produced 
turbulence and pressure, on the other hand. Here the amplitude 
of the primordial fields varies from B, = 0.05 to 0.50nG. The 
first two columns correspond to models with relatively weaker 
magnetism. The gas density maps in these two models resemble 
their counterparts in the first row. Indeed, the impact of weak pri- 
mordial fields on the matter power spectrum and consequently 
the buildup of galaxies is slim (see also Sethi & Subramanian 
2005; Sanati et al. 2020; Katz et al. 2021). However, the average 
gas mass experiences a slight increase, which can be attributed 
to the additional perturbations induced by primordial magnetic 
fields in the matter distribution. These perturbations lead to an 
earlier formation of the first dark matter minihalos and an earlier 
onset (approximately 100 Myrs) of gas collapse, which subse- 
quently continues at a higher rate. 


Increasing B; introduces considerable reformations in the 
main halo of the galaxy and the number of its subhalos. Dark 
matter only simulations of Sanati et al. (2020) shows a strong 
increase in the maximum of the Gaussian halo mass function, 
up to a factor of 7, when Ba = 0.20 and 0.50 nG, compared to 
their model without magnetically perturbed initial density fluctu- 
ations. For B, = 0.20 nG, Fig. 1 suggests that the additional per- 
turbations generated by primordial magnetic fields mainly affect 
halos with the Jeans mass M; x 108 Mo. This impact is reflected 
in Fig. 2, where the number of subhalos in this mass range is 
noticeably increased. Moreover, Mgas inside the galactic region, 
which contains the main halo of the targeted dwarf and its sub- 
halos, is increased. 
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Fig. 2: Gas density projections of model galaxy h070. The strength of the magnetic field is increasing from B; = 0.05 nG on the 
left to Ba = 0.50nG at the rightmost panels. The upper and lower rows represent all the models generated without primordial 
magnetic fields (ACDM), and with primordial magnetic fields (ACDM+PMFs) matter power spectrum, respectively. Snapshots are 
taken at redshift z = 9. The white circle shows the galactic region, which is traced in all time steps. The evolution of various galaxy 
quantities, such as the gas content, is evaluated inside this region with radius ~ 1.5 virial radius of the dark matter halo. Panels have 
physical sides and depth of 256 kpc. The included gas mass represents the average mass of gas in the galactic region over the first 


2 Gyr. 


When B; = 0.50nG, from Fig. 1 it appears that the altered 
matter power spectrum influences the formation of halos with 
mass Mj 10? Mo. Therefore, in this model, the main halo of 
the dwarf galaxy is primarily affected. In the final panel of Fig. 2 
it is evident that the gas mass is approximately 100 times greater 
compared to the counterpart derived from the ACDM paradigm 
depicted in the first row. Here, the strong gravitational potential 
of the central halo, can merge the satellite subhalos and leads to 
the formation of a massive dwarf galaxy. 


4.2. The evolution of energy components 


In this section, we explore the co-evolution of the key compo- 
nents shaping the energy allocation as a dwarf galaxy evolves 
over time. The three panels in Fig. 3, respectively from top to 
bottom, represent the time evolution of the magnetic energy, in- 
ternal energy and our estimate for the compressional evolution of 
the magnetic energy. For this plot we use dwarf spheroidal h878, 
but the outcomes are similar for other halos examined in this 
work (see Table 2). Note that the gas content within the galaxy 
halo, at different time steps, is subject to the natural mass growth 
due to accretions, mergers and outflows. Equally, computing the 
total energy over time is sensitive to the shape and volume of 
a galaxy. Therefore, for a fair comparison between different en- 
ergy components, we use specific energies 


E, 
Maas í 


(16) 


£y = 


as the ratio between the total energy Ex of component x and the 
gas mass Mgas, inside the galactic region (see Section 3.4). The 
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main amplification of the magnetic energy is expected to occur 
before z — 2 (Bernet et al. 2008). During this period, gas kinetic 
energy is converted into magnetic through magnetic field ampli- 
fication by a small-scale dynamo, reaching equipartition levels 
(Groves et al. 2003). After this redshift, magnetism is observed 
to be in equipartition with other energy components (Wolfe et al. 
1992; Basu & Roy 2013; Beck 20152). For that reason and for 
clarity, we truncate the time evolution of energies in this plot 
at 2 Gyr. Each line represents a different model from Table 1, 
where color codes are the same as in Fig. 1. In the left column, 
the galaxy simulated using a classical matter power spectrum 
with the strength of the magnetic field varying from B; = 0.05 
to 0.50 nG, showcases B05, B10, B20, and B50 models. In the 
right column, the magnetically-induced matter perturbations by 
primordial fields are also included in the models. 


Magnetic energy The topmost panels show the specific mag- 
netic energy Emag as it evolves over time. The time evolution of 
the magnetic field is described by the induction equation: 


lódB) _ 
ð 


1 
V x |v x (a° B1) - -17 x (à B)|. (17) 
a 


a 


In the absence of resistivity (7 = 0), peculiar velocities (v = 0) 
and non-ideal magnetic sources such as Biermann battery (Bier- 
mann 1950), this leads to 0 (a? B;)/Ót = 0. Before density fluc- 
tuations start collapsing into structures, both magnetic fields and 
gas density are decreasing with the expansion of the Universe as 
1/a? and 1/à?, respectively, shaping the relation B oc 5 The 
magnetic flux freezing to the plasma is imprinted in the early 
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Fig. 3: Comparison of different specific energy components within the galactic region (white circle in Fig. 2) of h878. The left and 
right panels represent all the MHD runs emerging from a standard ACDM and ACDM+PMFs matter power spectrum, respectively. The 
upper panels show the time evolution of specific magnetic energy. The ratio of magnetic energy between different models with initial 
seeds between 0.05 to 0.50 nG increases by a factor smaller than BZ /87, suggesting energy loss due to magnetic field saturation, 
including energy dissipation during the collapse. The middle panels show the specific internal energy, which tracks the onset and 
rate of star formation. The lower panels show the estimate for the compressional part of the magnetic energy from isotropic adiabatic 
gas collapse. Magnetic energy initially follows the compressional energy as B oc B until magnetic field amplification causes a 


deviation in their time evolution. 


decaying of the specific magnetic energy at high redshifts, as 
shown in the topmost panels of Fig. 3. 


As expected, models with stronger initial seeds (B; = 0.20 
and 0.50 nG) reach higher magnetic energy values during their 
evolution. This pattern is evident in both the unperturbed ACDM 
set of simulations and the magnetically-perturbed set, as de- 
picted in the left and right columns, respectively. 


It is interesting that the energy ratios in each set of runs do 
not increase by a factor of B? /87. For instance, when we com- 
pare the B0.05 and BO.10 models, one might anticipate a 2? 
increase in energy due to a stronger initial magnetic seed. How- 
ever, a smaller energy increase of ~ 2.7 suggests energy loss due 
to magnetic field saturation, with some energy dissipation even 
occurring during the collapse phase. The saturation of magnetic 
energy becomes more pronounced in stronger models and can 
occur within relatively short time intervals, such as a few hun- 
dred Myrs. This behaviour is followed in both sets of simulations 
presented in the left and right panels. In models featuring weaker 


magnetization (B; = 0.05 and 0.10 nG), it takes a longer period 
for the magnetic energy to reach saturation. 


Internal energy The middle panels show the specific internal 
energy £in, aS a function of time. In the ACDM setup, the in- 
ternal energy gradually decreases as the value of B; increases 
from 0.05 to 0.50 nG. This decline can be attributed to the sup- 
plementary pressure support exerted by magnetic fields, which 
also opposes the collapse of gas clouds (see the discussion in 
Section 4.1) and hinders star formation. Given that the internal 
energy closely tracks the heating generated by the stellar feed- 
back, raising the magnetic field strength leads to a decrease in 
the number of peaks and valleys, as well as the amplitude of int. 

This reduction in internal energy is particularly pronounced 
during the early stages of galaxy formation, specially at redshifts 
z > 6, distinguishing the stronger models B8.28 and BO.50 
from the weaker models B8.85 and BO. 18. In the extreme case 
of BQ.50, the star formation is dramatically suppressed and a 
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limited increase of £int is observed. Another notable charac- 
teristic evident here is the delay in the growth of internal en- 
ergy when the magnetic field strength is significant, particularly 
B; = 0.20 and 0.50 nG. The evolution of internal energy closely 
corresponds to changes in density and temperature and conse- 
quently reflects the collapse of gas clouds. Therefore, the delay 
caused by strong magnetic fields in the formation of minihalos 
impacts the timing at which £in begins to accumulate from its 
initial value. 


In the ACDM+PMFs setup, a distinct trend is evident in the right 
panel compared to the left panel, showcasing the influence of pri- 
mordial fields. Here, a higher variability is evident in the internal 
energy, characterized by numerous peaks and valleys, specially 
in stronger models. This pattern, visible across all strengths of 
magnetic field, implies a more dynamic star formation history. In 
these models, the formation of more massive minihalos acceler- 
ates the pace of collapse, star formation and consequently influ- 
ences the rate at which internal energy is growing. Additionally, 
primordial magnetic fields expedite the onset of minihalo forma- 
tion. As a result, the collapse of gas clouds is accelerated, leading 
to the onset of internal energy accumulation within timescales as 
short as several tens of Myr, specifically when B; = 0.50 nG. 


Compressional energy The bottom panels show the evolution 
of the specific compressional energy €comp versus time. This 
term is obtained from an isotropic adiabatic approximation of 
the gas collapse. During the early stages of galaxy formation, 
the magnetic energy density scales with the gas density as œ Pis 
(see e.g., Rieder & Teyssier 2017; Martin-Alvarez et al. 2018). 
Thus Ecomp(to) = Emag(to) œ pgs (to) X volume(to). In the ab- 
sence of amplification or decaying processes, this approxima- 
tion dictates the evolution of the frozen-in magnetic field lines, 
which closely follow the collapse of cosmic structures, indicat- 
ing Eeomp(t) « pap x volume(f). For extracting this term in 
each time step of the simulation, we use 


dx Pgas cell (£) = 
Ecomp,cell(t) = ae (Emaggal(to)) X (gaa , (18) 
° 37r, (fo) = (Pgas.gal(to)) 


where the compressional energy in each cell of size dXce and 
at each time step t, denoted as Evompcen(t), is directly propor- 
tional to the gas density pgascen(/) within that particular cell. 
(Emag,gal(to)) stands for the average of magnetic energy within 
the galactic region of radius rga; at time tọ when it is at its min- 
imum state. This corresponds to the turnaround point, marking 
the moment when the initial perturbations transition from ex- 
pansion to collapse. At this moment, when the magnetic energy 
reaches its lowest value, baryonic matter accumulates and the 
first stars begin to form within the emerging galaxy. Eventually, 
the total compressional term is derived by summing over all cells 
within the galactic region. It is then divided by the total gas mass 
in this region to obtain the specific energy. By comparing the 
evolution of magnetic energy and compressional energy, we can 
investigate the nature of magnetic field amplification. In case of 
dynamo amplification, it is anticipated that magnetic energy un- 
dergoes exponential growth, and the ratio of magnetic energy to 
compressional energy significantly exceeds 1. 

Figure 4 presents the ratio between Emag and Ecomp as a func- 
tion of time. As described, the compressional term demonstrates 
the expected buildup of magnetic fields, when exclusively frozen 
in the collapse of matter into dark matter halos. The deviation of 
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Emag from Ecos, apparent in all models, indicates the amplifica- 
tion of magnetic fields through different mechanisms that occur 
alongside adiabatic magnetic compression. These mechanisms 
include the activity of a dynamo within collapsing gas clouds. In 
both upper and lower panels, the ratio initially follows the same 
curve in all models, as magnetic energy evolves with the expan- 
sion of the Universe as By œ a?, until the density perturbations 
start to collapse. Once the galaxy reaches the turnaround point 
from expansion to collapse in each model, the ratio takes dif- 
ferent paths. Increasing the magnetic field strength results in a 
slower amplification rate. This is due to a considerably higher 
By(to) and B; in the galaxy. Thus it takes significantly longer 
time for the stronger models to reach a similar amplification rate 
as the weaker ones. The upper and lower panel showcase the 
same result in ACDM and ACDM+PMFs setting. 

It is worth noting that this amplification can also be at- 
tributed to the influence of powerful supernova feedback events. 
In this scenario, the supernova shock evacuates the surround- 
ings of the exploding star from gas, thus reducing the Ecomp. 
Simultaneously, the supernova winds push the magnetic fields 
outwards, while keeping them in the galactic region. Therefore 
the Emag/Ecomp ratio is increasing with time. 

In order to unravel the source of separation between the com- 
pressional and magnetic energy terms, we re-simulate h070 in 
model B05 and B20, deactivating the supernovae feedback. The 
comparison is shown in the upper panel of Fig. 4. The dotted 
lines represent the simulations in which the supernovae feed- 
back is set to zero. The more rapid increase in magnetic en- 
ergy compared to compressional energy remains evident in sim- 
ulations without feedback. This indicates that the amplification 
of magnetic fields is not solely attributed to the supernova ex- 
plosion process. It further strengthens the notion that magnetic 
fields initiated with even a week seed, e.g., By = 0.05 nG can 
undergo dynamo amplification as the galaxy evolves through 
gas accretion, halo collapse, and mergers. The higher ratio of 
Emag/Ecomp in models without supernovae feedback is because 
magnetic field back reactions from supernovae explosions slow 
down the growth of magnetic energy in models with feedback. 


4.3. Impact of magnetic fields on star formation 


In this section, we present the influence of primordial magnetic 
fields on the star formation history and final stellar mass of dwarf 
galaxies. From each zoom-in simulation, we first identify the 
halo corresponding to the reference halo in the fiducial model. 
Subsequently, using the position and particle IDs, we extract the 
correspondent halo in other models. The physical quantities of 
each galaxy such as the stellar mass are computed inside the 
virial radius R»oo of the extracted halo. 

Figure 5 displays the cumulative mass of stars forming 
within the initial 2 Gyr. In each magnetic field model, the line 
represents the average value of stellar mass measured for all 
galaxies simulated using the same magnetic field strength. In 
the fiducial model B86, shown in dark blue, star formation is 
quenched in the simulated galaxies before 1 Gyr at the latest. 

In the first set of simulations, the delay induced in the onset 
of star formation in the MHD configuration can be seen in the 
top panel of Fig. 5 across all amplitudes of magnetic field exam- 
ined. This effect is specifically evident in the first Gyr. Follow- 
ing the delay in the birth of first stars, it is difficult for galaxies to 
catch up with the cumulative star formation in the fiducial model, 
when the gas is heated by the background UV-radiation. This 
delay becomes more apparent in models with strong magnetic 
fields, B; = 0.20 and 0.50 nG, which have a more substantial 
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Fig. 4: Time evolution of the ratio of magnetic energy to com- 
pression energy in the galactic region (white circle in Fig. 2) for 
all the simulations during the first 3.5 Gyr. The upper and lower 
panels corresponds to all the models in ACDM and ACDM+PMFs 
settings, respectively. The solid and dotted lines correspond to 
the runs with the fiducial value of supernovae feedback Esn = 
10?! erg and zero feedback, respectively. 


influence on the gas distribution. In weaker models, B; = 0.05 
and 0.10 nG, this delay results in a shortened star formation his- 
tory of less than 500 Myr, particularly pronounced in ultra-faint 
dwarf galaxies with a small gas mass. The rate of star forma- 
tion and the final stellar mass are consistently decreasing as the 
strength of the magnetic field increases. 


The bottom panel of Fig. 5 shows models which include 
magnetic fields through the altered matter power spectrum 
ACDM-PMFs. One striking outcome of these models is an expe- 
dited onset of star formation. This stems from the influence of 
primordial fields in speeding up the formation of first dark mat- 
ter halos and their subsequent gas accretion. This effect persists 
across all amplitudes of magnetic fields examined, further con- 
firming the results obtained in Sanati et al. (2020), for different 
amplitudes and spectral indices of primordial magnetic fields. 
However, it is more moderate in models with a weak primor- 
dial field. From Fig. 1 it can be seen that the region affected by 
weakly perturbed models (BO. 95ps and BO. 10ps) encompasses 
halos with a Jeans mass of approximately 10? — 107 Mo. Conse- 
quently, the weakly perturbed models exhibit a slight increase in 
the abundance of small minihalos within this mass range. How- 
ever, the gas content in most of these minihalos is insufficient 
to form stars. In these models, the accelerated build up of the 
galaxies is therefore a secondary effect of primordial magnetic 
fields when considering stellar masses. Their dominant effect is 
further counteracting pressure in the ISM. 

Increasing the amplitude of primordial magnetic fields af- 
fects the matter power spectrum at higher halo masses, as shown 
in Fig. 1. When B; = 0.20 nG, the abundance of subhalos with a 
Jeans mass of approximately 10° Mo increases. The gas collapse 
in these halos is able to satisfy the MTT prescription for the for- 
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Fig. 5: Cumulative stellar mass during the first 2Gyr. Each 
dashed line represents the stellar mass of each model galaxy 
simulated in a magnetic field model with strength varying from 
B; = 0.05 to 0.50 nG, compared with the model without mag- 
netic field. Each solid line represents the mean stellar mass of 
all galaxies simulated in each model. The upper and lower pan- 
els correspond to all the models in ACDM and ACDM+PMFs set- 
tings, respectively. Top) Magnetic fields delay the formation 
of first stars. In low mass halos of ultra-faint dwarf galaxies, 
Ba = 0.50nG completely prevents the star formation. In order 
to include halos which do not form any stars in the plot, their 
stellar masses are shown at 4 - 10? Mo. Bottom) The influence 
of primordial magnetic fields on the matter power spectrum ac- 
celerates the formation of the first dark matter halos, leading to 
an earlier onset of star formation but an overall reduction in the 
final stellar mass. 


mation of stellar particles. The star formation rate, specially dur- 
ing the first 500 Myr, is therefore considerably boosted in this 
model. Similarly, when B, = 0.50nG, the main halo of each 
dwarf galaxy grows in mass. Unlike lower amplitudes of B;, for 
these two models the impact of primordial magnetic fields in ac- 
celerating galaxy evolution outweighs their contribution to the 
overall pressure support and ISM physics. Figure 5 shows that 
in these two strongly perturbed models star formation is initi- 
ated earlier and at higher redshifts compared to the fiducial BOO 
model. This implies that primordial magnetic fields can be a po- 
tential mechanism to explain the efficient and rapid star forma- 
tion observed in galaxies at z ~ 8.5 — 14.5 in the CEERS survey 
(e.g., Finkelstein et al. 2023; Asada et al. 2024). 


To quantitatively demonstrate the correlation between star for- 
mation and magnetic field strength, Fig. 6 presents the final stel- 
lar mass of dwarf galaxies as a function of By. In the left panel, 
each circle is representative for one galaxy at redshift z = 0, 


Article number, page 11 of 16 


A&A proofs: manuscript no. main 


Mago [M ə] 


e ACDM + ACDM + PMFs 
10°F 
= è 
= o 3 E e 
E e e $ + 10° 
= 10 € = e i T $ $ 
E e e ° 1 
z [e 3 + g T 
e. e. e * 
+ * + 
° e * 
e * 
10" 2: + e ACDM 
no stars e e e @ ACDM+PMFs * 
0.0 0.05 0.10 0.20 050 0.0 0.05 0.10 0.20 0.50 
By [nC] By [nC] 


Fig. 6: Stellar mass of each galaxy at redshift z = 0 for varying primordial magnetic field strength, color-coded by the virial mass 
of the dark matter halo. In the left panel, each circle represents one galaxy simulated in the MHD suite (ACDM model). In the 
right panel, galaxies shown as empty diamonds (ACDM+PMFs model) consider both magneto-hydrodynamics and the contribution 
of primordial fields to the matter power spectrum. The stellar mass of halos which do not form any stars is shown at M, < 10^ Mo 
instead of 0 Mo. Halos that do not reach redshift z = 0, due to extensive CPU-hours required , are identified by marking them with 


an upward arrow. 


color-coded by the virial mass of its dark matter halo. The gray 
line represents the average stellar mass for all simulated galax- 
ies as a function of magnetic field strength. For both ultra-faints 
and dwarf spheroidal, the final stellar mass remains largely unaf- 
fected in the weak models with B, = 0.05 nG. However, compar- 
ing the dwarf spheroidal to ultra faint dwarf galaxies in stronger 
magnetic field models shows that the low-mass gas content of 
ultra faints is very sensitive to any perturbations, namely, the im- 
pact of magnetic fields. Although the merging history and the 
star formation history in each ultra-faint galaxy is different, the 
average trend shows a steady decrease in the final stellar mass. 


In the BO. 20 model, only the dwarf spheroidal galaxy within 
our sample has a non-zero final stellar mass, while the collapse 
of minihalos and star formation in all of the simulated ultra- 
faint dwarfs is completely hindered by a magnetic field with a 
strength of B; = 0.20nG. For illustrating these halos in Fig. 6, 
their zero stellar mass is replaced by random values between 5 
and 6 x 10°Mo, and marked by no stars sign. A stronger mag- 
netic field with a strength of B; — 0.50 nG prevents the forma- 
tion of stellar particles with our MTT star formation model in 
all the ultra-faint dwarf galaxies. Here, the most massive dwarf 
spheroidal in our sample, h825, is used to obtain a threshold for 
the mass of dwarf galaxies where the impact of magnetic fields 
become more important in their evolution. Due to computational 
constraints, simulations for this galaxy are limited to the fiducial 
model B®.00 and the extreme MHD model B8. 58. By compar- 
ing this more massive dwarf with lower mass ultra-fain galaxies, 
it appears that magnetic fields severely quench the star forma- 
tion in dwarf galaxies with halo mass and stellar mass below 
Ma0 ~ 2.5 + 10? Mo and M, ~ 3- 106 Mo, respectively. 

The right panel of Fig. 6 shows models which include the im- 
pact of primordial magnetic fields on the matter power spectrum 
(ACDM+PMFs) in the MHD setup. Notwithstanding the modifica- 
tions occurring in the internal structure of dwarf galaxies at high 
redshifts, the impact of magnetic fields in the ISM remains a vi- 
tal factor in their evolution throughout the star-forming phase. 
The result of this dual effect on the final stellar mass of model 
galaxies is displayed in colored empty diamonds. In strongly per- 
turbed models, B20ps and B5Qps, despite an increased number 
of minihalos with a larger gas content, the persistent effect of ad- 
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ditional magnetic energy in the ISM counteracting gravitational 
collapse throughout the entire star formation period leads to a 
lower stellar mass compared to the fiducial BOO model at red- 
shift z = 0. The noteworthy observation here is that the final 
stellar mass in halos emerging from the second set of simula- 
tions (ACDM+PMFs) surpasses that of their counterparts in the 
first set (ACDM). Therefore, for dwarf galaxies with luminosity 
Ly < 10° Lo, solely incorporating magnetic fields in numerical 
simulations while disregarding the influence of primordial fields 
on their host dark matter halos may result in miscalculation of 
the impact of primordial magnetic fields in reducing the star for- 
mation rate, leading to an overestimation for stronger magnetic 
field strengths. 


4.4. Impact of magnetic fields on scaling relations 


In this section, we explore the impact of primordial magnetic 
fields on the global properties of dwarf galaxies. Namely, the 
V-band stellar luminosity, stellar mass, virial mass, line-of-sight 
velocity dispersion, half-light radius, and stellar metallicity. The 
method for calculating each quantity is described in Section 3.4. 


4.4.1. Stellar mass-halo mass 


Figure 7 shows the relation between galaxy stellar mass and 
dark matter halo mass for each simulated dwarf galaxy at red- 
shift z = 0. The results from the abundance matching models of 
Behroozi et al. (2013) and Moster et al. (2013) are represented 
by solid and dashed lines, respectively. Both models have their 
halo masses extrapolated to 10° Mc to encompass the mass range 
of the simulated ultra-faint dwarf galaxies. In the fiducial model 
BOO, the resulting stellar masses range from approximately 10^ 
to 2.5 - 10° Mo. All galaxies in this model lie well within the ex- 
pected values for the final halo mass, which spans from 8 - 108 
to 2.5 - 10? Mo. Note that the stellar masses of these galaxies 
are lower than that those of the galaxies simulated in Revaz & 
Jablonka (2018b) and Sanati et al. (2020), despite having the 
same halo mass. This is due to the star formation and stellar 
feedback schemes varying in different hydrodynamical codes. 
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Fig. 7: Relation between galaxy stellar mass and dark matter 
halo mass for each simulated dwarf galaxy at redshift z = 0. 
For comparison, the extrapolated abundance matching models 
of Behroozi et al. (2013) and Moster et al. (2013) are shown 
in solid and dashed lines, respectively. In order to include halos 
which do not form any stars in the plot, their stellar masses are 
shown at M, ~ 10? Mo. 


Colored circles show halos in the first set of simulations, 
originating from a ACDM paradigm. Magnetism of low strength 
(Bı = 0.05 and 0.10nG) does not dramatically influence their 
location in the M, — Moo relation compared to the fiducial 
model. Nonetheless, the inverse scaling of the stellar mass with 
the strength of the field, observed in Section 4.3, persists in this 
relation as well. For a given halo mass, the final stellar mass de- 
creases with increasing B,. Low mass halos which do not form 
stars in the strong magnetic field model, where Ba = 0.20 and 
0.50 nG, are labeled with no stars. To include these halos in the 
plot, their zero stellar mass is substituted with a random value of 
5—-6x10?Mgs. 


Colored empty diamonds depict halos in the second set 
of simulations, emerging from a modified ACDM paradigm 
(ACDM+PMFs). The global properties of galaxies undergo only 
minor changes in the weakly perturbed models. When B; = 0.05 
and 0.10 nG, all halos closely resemble their counterparts in the 
fiducial model. However, in the strongly perturbed models, the 
final stellar mass recurrently reveals that the influence of primor- 
dial magnetic fields on the initial phases of galaxy formation, 
primarily through modifications to the matter power spectrum, 
outweighs its subsequent effects in the ISM. Galaxies which do 
not form stars in the BO. 20 and BO. 56 model, evolve into bright 
dwarfs with average M, = 6 and 7 - 10^ Mo in the BO. 20ps and 
BO.50ps models, respectively. Since the abundance matching 
relation is not well constrained in the regime of ultra-faint dwarf 
galaxies, none of the models is ruled out based on this compar- 
ison. It is indeed likely that the same halo mass supports orders 
of magnitude differences in stellar mass, owing to various envi- 
ronmental effects and unique history of each dwarf galaxy (see 
ie. Revaz & Jablonka 2018b). This may lead to a considerable 
scatter in the low mass regime. 
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Fig. 8: V-band luminosity as a function of line-of-sight veloc- 
ity dispersion in each model galaxy is compared to the observa- 
tional data of Milky Way satellites in gray squares. In the ACDM 
setting, increasing the magnetic field strength results in fainter 
systems (illustrated by colored circles) and can even completely 
prevent star formation in ultra-faint dwarf galaxies. In order to 
include halos which do not form any stars in the plot, their lu- 
minosities are represented at Ly ~ 10? Lo. In the ACDM+PMFs 
setting, despite lower luminosities, velocity dispersion of galax- 
ies (illustrated by colored empty diamonds) increases due to a 
higher frequency of mergers. 


4.4.2. Luminosity-velocity dispersion 


Figure 8 shows the V-band luminosity Ly versus the line-of-sight 
velocity dispersion oLos for each simulated galaxy at redshift 
z = 0. The comparison to the observational data is obtained from 
the regularly updated Local Group and Nearby Dwarf Galaxies 
database of McConnachie (2012b), represented by gray squares 
in the figures. 


In the first set of simulations (ACDM), represented by col- 
ored circles, higher values of B; result in a reduced final stellar 
mass for a given halo mass. Therefore, increasing the strength 
of magnetic fields, causes the systems to become noticeably 
fainter than their corresponding analogues in the fiducial model, 
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Fig. 9: Half-light radius as a function of V-band luminosity of 
each model galaxy at redshift z = 0 is compared to Milky 
Way and M31 satellites in black squares, given by McConnachie 
(2012b). In order to include six halos which do not form any 
stars in the plot, their luminosities are shown at Ly ~ 10? Lo. 


or even completely stops the star formation. Additionally, mag- 
netic fields impact the line-of-sight velocity dispersion in each 
halo. In galaxies where magnetic fields of strength B; = 0.05 nG 
lead to less concentrated star formation and an increased half 
light radius (see Sec. 4.4.3), there is an associated increase of 
the velocity dispersion. 

In the second set of simulations (ACDM+PMFs), represented 
by colored empty diamonds, incorporating primordial magnetic 
fields through the matter power spectrum, both the luminosity 
and the stellar velocity dispersion of the simulated galaxies are 
modified. Particularly, in the strongly perturbed models, B20ps 
and B56ps, despite lower luminosities, the line-of-sight veloc- 
ity dispersion of each halo is notably larger than in the fiducial 
model. In these models, due to the impact of primordial mag- 
netic fields on the evolution of dark matter substructures, stars in 
the main halo of each galaxy migrate from numerous minihalos. 
In addition, because of greater number of minihalos forming in 
these models, there is increased frequency of mergers between 
minihalos. This results in the observed differences in the dynam- 
ics of the simulated galaxies. 


4.4.3. Size-Luminosity 


Figure 9 displays the half-light radius r;;? of the model dwarf 
galaxies as a function of their V-band luminosity at redshift 
z = 0. The Local Group dwarfs compiled from McConnachie 
(2012b) are represented by gray squares with error bars. The 
galaxies simulated in the fiducial model cover a wide range in 
luminosity and size. The half-light sizes in this model span ap- 
proximately from 300 to 850 pc. 

Galaxies in the first set of simulations, ACDM, are represented 
by colored circles. Here, the influence of magnetic fields corre- 
lates with the half-light radius of the simulated galaxies. The 
more compact the galaxies, the less effect magnetic fields have 
on their gas distribution and, consequently, on their star forma- 
tion. A magnetic field strength of B, = 0.05 nG in the majority 
of galaxies in our sample increases the half-light radius by an 
average value of 25 percent. However, for B; > 0.05 nG, mag- 
netic fields hinder star formation in small minihalos and by pre- 
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venting the formation of stars in the outskirts, shrink the size 
of ultra-faint dwarf galaxies. It demonstrates that incorporating 
magnetic fields of strength B; > 0.05 nG and B; « 0.50nG in 
hydrodynamical simulations leads to a reduction in the size of 
dwarf galaxies, offering a plausible mechanism for the formation 
of very compact ultra-faint dwarfs (See e.g., Revaz (2023) for 
a detailed discussion on reproducing compact ultra-faint dwarf 
galaxies). 

Simulated dwarfs in the second set of simulations, 
ACDM+PMFs, are represented by colored empty diamonds. In 
these models, the stellar half-mass is more widely dispersed 
within the halo of each galaxy. This is a direct consequence of 
the primordial magnetic fields, which increase the stellar line- 
of-sight velocity dispersion. Here, stars form in more abundant 
subhalos and will be brought to the main halo of the galaxy 
by mergers. This results in an increased kinematic energy com- 
ponent and higher velocity dispersion (Revaz 2023). The large 
angular momenta of these merging subhalos imply larger sizes 
for the simulated dwarf galaxies. This effect becomes especially 
pronounced when B; = 0.20 and 0.50 nG. 


To summarize, the z = 0 properties of halos formed in our 
fiducial model, provide a good representation of spheroidal and 
ultra-faint dwarf galaxies in the Local Group, serving as a ro- 
bust basis for predicting the influence of magnetic fields. From 
our findings, we can affirm that primordial magnetic fields with 
strengths B, > 0.05 nG significantly impact the physical observ- 
ables of dwarf galaxies, specially the small mass ultra-faints with 
luminosities below 10? Lo. In this work, the best consistency 
with the scaling relations of Local Group dwarfs in our simu- 
lated galaxies is obtained when incorporating primordial mag- 
netic fields with moderate strengths of B; = 0.05 and 0.10 nG in 
the MHD configuration. 


5. Discussion & Conclusions 


In this study, we investigated the impact of primordial mag- 
netic fields on the formation and evolution of dwarf galaxies. We 
conducted new cosmological zoom-in magneto-hydrodynamical 
simulations from redshift z — 200 to zero. These simulations are 
focused on eight halos that host dwarfs with V-band luminosities 
ranging from approximately 10? to 10° Lo. All simulations are 
generated using our modified version of the constrained trans- 
port MHD RAMSES code. 

We employed two different simulation setups. In the first 
setup, a uniform magnetic field is initialized with a given 
strength and the matter power spectrum follows the classical 
ACDM paradigm. In the second setup, we explored the impact of 
both varying magnetic field strengths and including the effect of 
a Gaussian random primordial magnetic field on the matter dis- 
tribution. This approach allowed us to examine the combined ef- 
fect of primordial magnetic fields on the matter power spectrum, 
on the one hand, and the influence of magneto-hydrodynamics 
on the evolution of small mass galaxies, on the other hand. 

Depending on the strength, each magnetic field model af- 
fects a distinct mass range of dark matter halos, within the mass 
range relevant to host halos of dwarf galaxies. In this study, 
we explored a range of primordial magnetic fields with varying 
strengths, specifically in the range B; = 0.05 — 0.50 nG, while 
keeping the spectral index constant at ng = —2.9. All of these 
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magnetic field strengths fall within the limits allowed by cosmo- 
logical constraints listed in Section 1. Our main results can be 
summarized as follows: 


— In the context of a ACDM cosmological framework, mag- 
netic fields with an initial strength of B; > 0.05 nG, have 
a non-negligible contribution in shaping the evolution of 
dwarf galaxies. These fields not only delay the collapse of 
gas clouds, but in the presence of Ba > 0.20nG prevent 
the formation of dense gas minihalos altogether. This ef- 
fect is stronger in low mass ultra-faint dwarf galaxies with 
halo mass and stellar mass below Mago ~ 2.5 - 10? Mo and 
M, ~ 3- 10° Mo, respectively, due to the shallow potential 
well of their host halos. 

— [n a cosmological setup evolved from initial conditions of 
ACDM+PMFs, the main halo of a dwarf galaxy, as well as 
the number of its subhalos, are influenced. When the halos 
affected by primordial magnetic fields fall within the Jeans 
mass range M; < 10? Me, the gas is distributed among a 
higher number of these low-mass subhalos. On the other 
hand, if the halo mass range affected contains the Jeans mass 
of the main halo of a dwarf galaxy (M; ~ 10? Mo), the gas 
content in the central halo increases. This strong gravita- 
tional potential can lead to the merging of satellite subhalos 
and ultimately increase the final mass of the host galaxy. 

— Following the delay caused by magnetic fields during the for- 
mation of first stars, the gravitational field of low-mass dark 
matter halos is unlikely to retain gas when it is pushed out by 
UV-radiation during the epoch of reionisation. This results in 
the depletion of the gas reservoir in minihalos and an earlier 
quenching of star formation. In our sample, the density and 
temperature of gas in all six ultra-faint halos remain below 
the MTT criteria to form stellar particles. However, the in- 
fluence of primordial magnetic fields on the matter power 
spectrum accelerates the formation of the first dark matter 
minihalos, leading to an earlier onset and a higher rate of 
star formation at redshifts z > 12. 

— The observable properties of dwarf galaxies at redshift z = 0 
are indicative of the influence of magnetic fields on their star 
formation history and kinematics. Comparing the two sim- 
ulation setups, we observe that dwarf galaxies forming in a 
ACDM+PMFs Universe with a moderate primordial magnetic 
field of Ba = 0.05 and 0.10 nG are more consistent with 
the scaling relations of Local Group observations. However, 
stronger magnetic fields lead to dwarf galaxies with larger 
size and higher velocity dispersion compared to the observed 
Milky Way satellites. 

— Increasing the strength of magnetic fields does not result in 
a proportional B^ 87 increase in magnetic energy. This sug- 
gests energy loss due to magnetic field saturation, with some 
energy dissipation even occurring during the collapse phase. 
The saturation of magnetic energy occurs more rapidly when 
B; = 0.20nG, within relatively short time intervals of a few 
hundred Myrs. When B; x 0.10nG, it takes a longer time 
for the magnetic energy to reach saturation. 

— In our simulations, the buildup of magnetic energy in time 
differs from the evolution of magnetic fields in an adiabatic 
spherical collapse with magnetic flux freezing, also referred 
to as compressional amplification. This latter is not driven by 
the supernovae feedback, suggesting that the amplification of 
the magnetic energy is not provoked in the course of super- 
novae explosions. Instead, it supports the growth of magnetic 
fields initialized with a seed as small as B; = 0.05 nG due to 
dynamo amplification. 


In summary, primordial magnetic fields have a substantial 
impact on the initial stages of dwarf galaxy formation by modi- 
fying the matter power spectrum. Moreover, magnetic fields play 
an important role in the evolution of dwarf galaxies by contribut- 
ing to the ISM energy budget. These effects situate dwarf galax- 
ies as an unprecedented window into primordial magnetism. 
Therefore, the interplay between simulations and observations of 
spheroidal and ultra-faint dwarf galaxy populations has the po- 
tential to further constrain and inform our understanding of the 
primordial magnetic fields in our Universe. On the other hand, 
the early structure formation induced by primordial magnetic 
fields is potentially capable of producing rapid and efficient star 
formation at high redshifts suggested by the JWST data. 
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